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To describe the microscopic properties of matter, quantum mechanics uses wave functions, whose 
structure and time dependence is governed by the Schrödinger equation. In atoms the charge distributions 
described by the wave function are rarely observed. The hydrogen atom is unique, since it only has one 
electron and, in a de electric field, the Stark Hamiltonian is exactly separable in terms of parabolic 
coordinates (n, é, œ). As a result, the microscopic wave function along the € coordinate that exists in the 
vicinity of the atom, and the projection of the continuum wave function measured at a macroscopic 
distance, share the same nodal structure. In this Letter, we report photoionization microscopy experiments 
where this nodal structure is directly observed. The experiments provide a validation of theoretical 
predictions that have been made over the last three decades. 
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The development of quantum mechanics in the early 
part of the last century has had a profound influence on 
the way that scientists understand the world. Central 
to quantum mechanics is the concept of a wave function 
that satisfies the time-dependent Schrödinger equation 
[1]. According to the Copenhagen interpretation, the 
wave function describes the probability of observing the 
outcome of measurements on a quantum mechanical 
system, such as measurements of the energy or the position 
or momenta of constituents [2]. The Copenhagen interpre- 
tation thus allows reconciling the occurrence of nonclass- 
ical phenomena on the nanoscale with manifestations and 
observations made on the macroscale, which correspond to 
viewing one of a number of possible realizations allowed 
for by the wave function. 

Despite the overwhelming impact on modern electronics 
and photonics, understanding quantum mechanics and the 
many possibilities that it describes continues to be intellec- 
tually challenging, and has motivated numerous experi- 
ments that illustrate the intriguing predictions contained 
in the theory [3]. Using ultrafast lasers, Rydberg wave 
packet experiments have been performed illustrating how 
coherent superpositions of quantum mechanical stationary 
states describe electrons that move on periodic orbits 
around nuclei [4]. The wave function of each of these 
electronic stationary states is a standing wave, with a 
nodal pattern that reflects the quantum numbers of the 
state. Mapping of atomic and molecular momentum wave 
functions has been extensively explored by means of 
(e, 2e) spectroscopy, using coincident detection of the 
momentum of both an ejected and a scattered electron to 
retrieve the momentum distribution of the former prior to 
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ionization [5]. In the spirit of scanning tunneling methods, 
orbital tomography based on high harmonic generation 
was developed as a method allowing the determination of 
atomic and molecular orbitals [6,7]. In this Letter we will 
present experiments where the nodal structure of electronic 
wave functions of hydrogen atoms is measured, making 
use of a photoionization microscopy experiment, where 
photoelectrons resulting from ionization after excitation 
of a quasibound Stark state are measured on a two- 
dimensional detector. 

The hydrogen is a unique atom, since it only has one 
electron and, in a dc electric field, the Stark Hamiltonian is 
exactly separable in terms of parabolic coordinates. For 
this reason, an experimental method was proposed about 
thirty years ago, when it was suggested that experiments 
ought to be performed projecting low-energy photoelec- 
trons resulting from the ionization of hydrogen atoms onto 
a position-sensitive two-dimensional detector placed per- 
pendicularly to the static electric field, thereby allowing the 
experimental measurement of interference patterns directly 
reflecting the nodal structure of the quasibound atomic 
wave function [8—10]. 

In a static electric field F the wave function of atomic 
hydrogen can be separated in terms of the parabolic 
coordinates 7, é, ¢ (n =r— z and é= r + z, where r 
is the distance of the electron from the proton, z is the 
displacement along the electric field axis and g = 
tan™!(y/x) is the azimuthal angle [see Fig. 1(a) and 
Ref. [11]]). Note that atomic units are used, unless speci- 
fied otherwise. Consequently, the wave function may be 
written as a product of functions y,(€) and y(n) that 
separately describe the dependence along é and 7, 
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FIG. 1. (a) Potential energy landscape and relevant coordinate 


system for hydrogen atom photoionization microscopy in an 
808 V/cm electric field (a.u. = atomic units). The hydrogen 
atom sits at the origin of the (z, y) coordinate system and the 
electric field is along the z axis. The boundary between the 
classically allowed and the classically forbidden region is plotted 
(solid lines) at the excitation energies of the four measurements 
that are shown in Fig. 3, i.e., ranging from E = — 172.8 cm™! to 
E = —163.3 cm™! (thick outer solid line). Close to the saddle 
point, the electron can only escape through a very narrow gap in 
the Coulomb + dc field potential. The parabolic coordinates 
n =r — z and é = r + z are illustrated by plotting a series of 
contours at constant 7 (dashed lines) and £ (dotted lines). The 
electron motion is always bound in the € coordinate whereas the 
motion along the 7 coordinate depends on the energy available 


for the 7 motion; (b) and (c) Potential ener y curves 
V(n)=- 2 ae 2—1 =i and V(€) = ae + a + F£ describ- 


ing the a ia the ņ and é coordinates Cit] where 
Z, =(n, + m/n and Z, = (n, + mt yy, V(m) and V(é) 
are shown for the (nı, n2, m) = (3, 26,0) quasibound state at 
E = —163.3 cm™!. 


ie, WE, n, o) = Qané) 7x1 lE x(n)e"?. The func- 
tions yı(é) and y(n) satisfy the ordinary differential 
equations, 
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In these expressions m is the magnetic quantum number 
and Z, and Z, are separation constants subject to the 
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condition Z; + Z) = 1. The parabolic quantum numbers 
nı and n, are related to the principal quantum number via 
n =n; +m + |m| + 1. Functions y,(€) and y(7) have 
nı and n, nodes along the é and 7 coordinates, respec- 
tively. Figure 1(a) shows the potential energy landscape 
for the hydrogen atom in an 808 V/cm electric field. The 
electron motion is always bound in the é coordinate 
whereas the motion along the 7 coordinate depends 
on the energy available for its motion. Figures in the 
lower panels illustrate the potential energy curves 
Vín) aN 2 + _ mn and V(é) = -3 + ve + Fe 
(please note that in ie literature other ka for 
the definition of V(7) and V(é) are in use as well; see, 
e.g., Ref. [12]), defining the motion along the 7 and é 
coordinates, respectively. Reference [8] contained a 
remarkable prediction for the special case where the 
atomic hydrogen photoionization involves the excitation 
of quasibound Stark states. In this case, where both n, and 
n, are good quantum numbers and the electron tunnels 
through the barrier in the potential energy curve associated 
with the 7 coordinate, the measurements should show a 
total of n, dark fringes, directly revealing an important 
signature of the Stark state involved. However, to date, this 
experiment was never performed. 

Motivated by the theoretical predictions for the con- 
figuration of the above-mentioned “photoionization 
microscope” [8], a photodetachment microscope for 
negative ions was first constructed by Blondel et al. [13]. 
Their experiments clearly revealed interferences between 
the photoelectrons en route to the detector, in agreement 
with simple semiclassical considerations by Du [14]. 
In photodetachment, the photoelectrons follow one of 
two possible parabolic trajectories to the two-dimensional 
detector. By contrast, in a photoionization experiment that 
starts from a neutral sample, the photoelectrons move in a 
combined static electric + Coulomb field, significantly 
complicating the dynamics and leading to the existence 
of an infinite number of trajectories that the photoelectron 
can follow to the detector [15]. 

Given the considerable challenges connected to the 
experimental use of atomic hydrogen, first attempts to 
implement photoionization microscopy were performed 
on Xe atoms by Nicole et al. [16]. Observed interference 
rings were interpreted in the framework of a semiclassical 
treatment [17], excluding the possibility of resonant exci- 
tation of a Stark state. The experiments themselves were 
performed both with and without the resonant excitation of 
quasibound Stark states, and no significant differences 
were observed [18]. More recently, photoionization mi- 
croscopy experiments were performed for Li atoms [19], 
revealing first indications of differences in the radial dis- 
tributions for on- and off-resonance excitation. This work 
provides the motivation for the hydrogen experiments 
reported in the present Letter, where we present results 
for both resonant and nonresonant ionization, and where 
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we convincingly validate the long time predictions of 
Demkov et al. [8,9]. 

An atomic hydrogen beam was formed by collimating 
hydrogen atoms resulting from the photodissociation of 
H,S gas in a first vacuum chamber (see Fig. 2 and the 
Supplemental Material [20]). The hydrogen atoms were 
ionized in the active region of a velocity map imaging 
(VMI spectrometer [21], where an 808 V/cm static elec- 
tric field was applied. The atoms were resonantly excited 
to a mixture of n = 2 s and p states by a two-photon 
transition (Aaser = 243 nm) and were ionized using 
narrowband, tunable laser pulses (Aaser = 365-367 nm, 
Tlaser = 8 ns) from a Fourier-limited, home-built pulsed 
dye amplifier [22]. The polarization of the 365-367 nm 
laser was along the static electric field (i.e., perpendicular 
to the detector). A dual microchannel plate (MCP) detector 
followed by a phosphor screen and a CCD camera was 
used to record the photoelectrons. An electrostatic zoom 
lens magnified the images by about 1 order of magnitude 
compared to the size that would have been measured 
without this lens [23]. 

The main results of the experiments are shown in Fig. 3. 
In this figure calculated and experimental results are shown 
for four experiments, where the hydrogen atoms were ex- 
cited to the (n,, n2, m) = (0, 29, 0), (1, 28, 0), (2, 27, 0) and 
(3, 26, 0) quasibound Stark states. As indicated in Fig. 3, the 
states lie at energies of — 172.82 cm™!, —169.67 cm“, 
—166.45 cm™!, and —163.30 cm™! with respect to the 
field-free ionization limit, i.e., just above the saddle point 
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FIG. 2 (color online). Schematic overview of the experiment. 
An atomic hydrogen beam was formed by photodissociating H3S 
and placing a 3 mm aperture (a) 65 mm downstream. In the 
active region of a velocity map imaging (VMI) spectrometer, the 
ground state hydrogen atoms were first excited to a mixture of 
n = 2 s and p states by a two-photon transition using a pulsed 
243 nm laser. Next, they were ionized by a Fourier-limited, 
tunable (365-367 nm), UV laser. By applying a voltage differ- 
ence across the repeller (b) and extractor (c) electrodes, the 
photoelectrons were accelerated towards a two-dimensional 
detector (d), consisting of a set of microchannel plates 
(MCPs), a phosphor screen and a CCD camera. En route to 
the MCP detector, the photoelectrons passed through a three- 
element Einzel lens (e), allowing an increase of the diameter of 
the recorded image by about one order of magnitude. 
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in the Coulomb + dc field potential, which lies at 
—174.00 cm™!. According to Eq. (6) of Ref. [24], the 
validity of which was checked experimentally [25], 
the ionization rate of these states covered a range from 
P=22x10!s! to T =7.25X 10°s™!, which 
[using SE(cm~!) = 5.3 X 107! T(s7!)] implies line 
widths comparable to the 0.005 cm! bandwidth of our 
excitation laser. These states could readily be identified 
in wavelength scans, since the ionization is complete 
before the hydrogen atoms leave the interaction region of 
the VMI [26]. By contrast, in the same energy range Stark 
states in the n = 31 manifold (T > 107 s~!) lead to very 
broad resonances, while states in the n = 29 manifold 
(T < 10° s~') undergo insufficient ionization before the 
atoms fly out of the interaction region. Total ionization 
spectra as a function of excitation energy in the given static 
electric field were successfully reproduced by means of the 
semiclassical Stark theory of Harmin [27,28]. The parabolic 
quantum number n, was identified by comparing the experi- 
mental spectra with separate theoretical excitation curves for 
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FIG. 3 (color online). Experimental observation of the 
transverse nodal structure of four atomic hydrogen Stark 
states. The images in the middle show experimental measure- 
ments for (nı, n, m) = (0,29,0), (1,28,0), (2,27,0), and 
(3, 26, 0). Interference patterns are clearly observed where the 
number of nodes corresponds to the value of nı. The results may 
be compared to TDSE calculations shown to the left (for details 
see text), revealing that the experimentally observed nodal struc- 
tures originate from the transverse nodal structure of the initial 
state that is formed upon laser excitation. A comparison of the 
experimentally measured (solid lines) and calculated radial 
(dashed lines) probability distributions P(R) is shown to the right 
of the experimental results. In order to make this comparison, the 
computational results were scaled to the macroscopic dimensions 
of the experiment. Please note that, since P(R) = f P(R, a)Rda, 
the radial probability distributions P(R) have a zero at R = 0, 
even if the two-dimensional images P(R, œ) do not. 
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each nı channel. Given the value of ,, the value of n was 
subsequently determined by applying the Wentzel-Kramers- 
Brillouin quantization rule along the 7 coordinate. 

The main result of the experiment, which is directly 
visible in the four images shown in the middle of Fig. 3, 
is the observation of an interference pattern with a number 
of dark fringes corresponding to the value of n,. This 
observation validates the prediction by Demkov and co- 
workers [8,9] and illustrates that photoionization micros- 
copy can be used to visualize the nodal structure of y;(€) 
for quasibound Stark states of the hydrogen atom. A ration- 
alization for this behavior can be found in the calculations 
shown to the left of the experimental images. Here, results 
of propagating the time-dependent Schrédinger equation 
(TDSE; see the Supplemental Material [20]) following 
excitation of the hydrogen atom at the energies used in 
the experiments are shown [29]. The graphs represent time- 
integrated plots of the two-dimensional electron density 
p(r, z) = r|V(z, t)|?, evaluated from the time of excitation 
(t = 0) to a time delay of 600 ps. The nodal structure that is 
observable at a large distance from the proton (here, 
0.4 um) clearly has its origin in the transverse nodal 
structure of the initial state that is formed upon laser 
excitation. We note that beyond a distance of 0.4 «wm the 
calculation shows no significant changes. 

A direct comparison of the experimental (solid line) and 
calculated (dashed line) results that is obtained by scaling 
the radial coordinate in the calculation, is shown to the 
right of the experimental results. Here a comparison of 
the measured radial probability distribution Poxp(R) = 
J P(R, œ)Rda [where P(R, a) represents the intensity 
distribution in the experimental image in polar coordinates 
R and a], and the calculated radial probability distribution 
Peatc(R) = R\Vcac(R, t)|? is given, showing very satisfac- 
tory agreement and validating the assignment of the num- 
ber of dark fringes as the parabolic quantum n,. Observed 
differences between the experimental and calculated 
results may be due to differences in the experimental 
conditions and the assumptions made in the calculations 
(where the Stark states were excited using a 250 ps exci- 
tation pulse), imperfections in the experimental images, 
and possible smearing effects due to the finite lifetime of 
the excited Stark states. 

A striking observation in the experiments is the pro- 
nounced difference between images recorded following 
resonant excitation of a quasibound Stark state and images 
recorded following nonresonant excitation of the ionization 
continuum. This is illustrated in Fig. 4, where a comparison 
is shown between the image for the (n,, no, m) = (2, 27, 0) 
state and two nonresonant images recorded 1.8 cm~! below 
and 1.1 cm~! above this resonance. Remarkably, the out- 
ermost ring in the on-resonance image extends significantly 
further radially than in the two surrounding images. In fact, 
when comparing the experimental results (triangles in the 
inset) to calculations using the classical method used in 
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FIG. 4 (color online). Evidence for on-resonance ionization by 
tunneling through the Coulomb + static field potential. A 
comparison is shown between a measurement carried out for 
the (nı, n,m) = (2, 27, O) resonance (b) and two non-resonant 
measurements performed 1.8 cm7! below (a) and 1.1 cm™! 
above (c) this resonance. The normalized radial distribution of 
the on- resonance measurement containing three maxima ex- 
tends significantly further outward than the two off-resonance 
measurements which show only a single maximum (d). The inset 
in (d) shows a comparison of the radial extension of the experi- 
mental images, defined as the position of the outer maximum 
(color triangles) and the theoretical radial extension (blue, solid 
line) according to the classical formula (excluding tunneling 
contributions) Rmax = [2L ajy 2 [15], where L is the distance 
between the H atom and the detector (L = 0.5 m), and E,, is the 
saddle point energy (— 174.00 cm™!). The experimental and 
theoretical radial extensions were matched for the measurement 
at E = —165.37 cm}. 


Ref. [15] (solid blue line in the inset), one sees that the 
position of the outer ring in the image for the (2, 27, 0) 
resonance extends further outwards by about 70%. This is 
in line with recent theoretical work by Zhao and Delos, who 
developed both a semiclassical and a quantum-mechanical 
theory for the hydrogen atom photoionization microscopy 
problem [30,31]. They predicted a “remarkable tunneling 
effect” that applies only in the case of resonant excitation 
of quasibound Stark states. Classically, the electron is 
trapped by the potential barrier V(7) [see Fig. 1(b)] if the 
emission angle is smaller than a critical angle 0, = 
arccos(1 — E?/2F). However, in case of excitation to a 
quasibound state, electrons with an emission angle smaller 
than 6. may tunnel through the V(7) potential barrier, 
leading to a situation where the electron can reach a posi- 
tion on the detector that is not classically accessible. 
Generally, the image measured at a resonance corre- 
sponds to a coherent superposition of resonant and 
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nonresonant contributions, the latter corresponding to 
direct excitation into the ionization continuum. As a con- 
sequence, a beating between these two contributions is 
expected. In the hydrogen measurements that are presented 
here, the resonant contribution strongly dominates. For 
example, the signal (i.e., the total number of detected 
electrons per acquisition) at the (2, 27, 0) resonance in 
Fig. 4 [case (b)] was stronger than the signal at the adjacent 
nonresonant positions [cases (a) and (c)] by a factor 10. 
Therefore, the image essentially represents a direct macro- 
scopic projection of the microscopic electronic quantum 
state. In other atoms electron-electron interactions 
(as manifested by quantum defects) that couple the initial 
state (nı, n2, m) to other states have a major influence on 
the electronic wave that is observed. For example, the 
above-mentioning tunneling in the 7 coordinate is largely 
absent in nonhydrogenic systems, because (nı, n, m) 
initial states that require tunneling couple to states that 
do not. In Xe, the coupling among parabolic states led to a 
complete disappearance of the resonant effects [16], 
whereas in Li it led to a substantial reduction of the 
contrast between resonant and nonresonant excitation [19]. 

In conclusion, we have shown that the concept of 
photoionization microscopy, as theoretically proposed 
more than 30 years ago and the subject of recent theoretical 
work predicting the possibility to experimentally observe 
nonclassical photoionization dynamics involving tunneling 
through the V(7) potential barrier, can be experimentally 
realized, providing both a beautiful demonstration of 
the intricacies of quantum mechanics and a fruitful play- 
ground, where the fundamental implications of this theory 
can be further explored. For example, predictions have 
already been made for the case where both electric 
and magnetic fields are present [32]. The experimental 
observations of the nodal structures of the wave functions 
presented in this Letter are not specific to the field strengths 
and laser excitation conditions (i.e., the polarization direc- 
tion of the exciting laser) used, but are a general phenome- 
non that is observable and can be exploited over a wide 
range of experimental conditions. 
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